Coherence and control of quantum registers based on electronic spin in a nuclear spin 

bath 



o 
o 

(N 

in 



43 

Oh 1 
i 

cr 



> 

q 

o 
o 



P. Cappcllaro 1,2 , L. Jiang 2 , J. S. Hodges 2 ' 3 and M. D. Lukin 1 ' 2 
1 IT AMP - Harvard-Smithsonian Center for Astrophysics and 2 Physics Department, 
Harvard University, Cambridge, MA 02138, USA. 3 Department of Nuclear Science and Engineering, 
Massachusetts Institute of Technology, Cambridge, MA 02139, USA 

We consider a protocol for the control of few-qubit registers comprising one electronic spin em- 
bedded in a nuclear spin bath. We show how to isolate a few proximal nuclear spins from the rest 
of the environment and use them as building blocks for a potentially scalable quantum informa- 
tion processor. We describe how coherent control techniques based on magnetic resonance methods 
can be adapted to these electronic-nuclear solid state spin systems, to provide not only efficient, 
high fidelity manipulation of the registers, but also decoupling from the spin bath. As an exam- 
ple, we analyze feasible performances and practical limitations in a realistic setting associated with 
nitrogen- vacancy centers in diamond. 



The coherence properties of isolated electronic spins in 
solid-state materials are frequently determined by their 
interactions with large ensembles of lattice nuclear spins 
[HQ- The dynamics of nuclear spins is typically slow, re- 
sulting in very long correlation times of the environment. 
Indeed, nuclear spins represent one of the most isolated 
systems available in nature. This allows, for instance, 
to decouple electronic spin qubits from nuclear spins via 
spin echo techniques 0, Q] • Even more remarkably, con- 
trolled manipulation of the coupled electron-nuclear sys- 
tem allows one to take advantage of the nuclear spin 
environment and use it as a long-lived quantum mem- 
ory 0, @, 0] • Recently, this approach has been used to 
demonstrate a single nuclear qubit memory with coher- 
ence times well exceeding tens of milliseconds in room 
temperature diamond Q. Entangled states composed 
of one electronic spin and two nuclear spin qubits have 
been probed in such a system The essence of these 
experiments is to gain complete control over several nu- 
clei by isolating them from the rest of the nuclear spin 
bath. Universal control of systems comprising a single 
electronic spin coupled to many nuclear spins has not 
been demonstrated yet and could enable developing of 
robust quantum nodes to build larger scale quantum in- 
formation architectures. 

In this Letter, we describe a technique to achieve co- 
herent universal control of a portion of the nuclear spin 
environment. In particular, we show how several of these 
nuclear spins can be used, together with an electronic 
spin, to build robust multi-qubit quantum registers. Our 
approach is based on quantum control techniques asso- 
ciated with magnetic resonance manipulation. However, 
there exists an essential difference between the proposed 
system and other previously studied small quantum pro- 
cessors, such as NMR molecules. Here the boundary be- 
tween the qubits in the system and the bath spins is ul- 
timately dictated by our ability to effectively control the 
qubits. The interactions governing the couplings of the 
electronic spin to the nuclear qubit and bath spins are of 
the same nature, so that control schemes must preserve 



the desired interactions among qubits while attempting 
to remove the couplings to the environment. The chal- 
lenges to overcome are then to resolve individual energy 
levels for qubit addressability and control, while avoiding 
fast dephasing due to uncontrolled portion of the bath. 

Before proceeding we note that various proposals for 
integrating these small quantum registers into a larger 
system capable of fault tolerant quantum computation or 
communication have been explored theoretically 1^, 11 1 
and experimentally 13, ■ These schemes generally re- 
quire a communication qubit and a few ancillary qubits 
per register, in a hybrid architecture. The communica- 
tion qubits couple efficiently to external degrees of free- 
dom (for initialization, measurement and entanglement 
distribution), leading to an easy control but also fast de- 
phasing. The ancillary qubits instead are more isolated 
and can act as memory and ancillas in error correction 
protocols. While our analysis is quite general in that it 
applies to a variety of physical systems, such as quantum 
dots in carbon nanotubes or spin impurities in silicon 
0, as a specific example we will focus on the nitrogen- 
vacancy (NV) centers in diamond 0, H, Q. These are 
promising systems for the realization of hybrid quantum 
networks due to their long spin coherence times and good 
optical transitions that can be used for remote coupling 
between registers [l5j |. 

To be specific, we focus on an electronic spin triplet, 
as it is the case for the NV centers. Quantum informa- 
tion is encoded in the to s ={0,1} Zeeman sublevels, sep- 
arated by a large zero field splitting (making m s a good 
quantum number). Other Zeeman sublevels arc shifted 
off- resonance by an external magnetic field B z , applied 
along the N-to-V axis. The electron- nuclear spin Hamil- 
tonian in the electronic rotating frame is 



h = lu l J2H + s z j: a 3 ■ P + H di] 



1-5, 



where S and P are the electronic and nuclear spin oper- 
ators, lot, is the nuclear Larmor frequency, Aj the hyper- 
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FIG. 1: a) System model, showing the electronic spin in red 
and the surrounding nuclear spins. The closest nuclear spins 
are used as qubits. Of the bath spins, only the spins outside 
the frozen core evolve due to dipolar interaction, causing de- 
coherence, b) Frequency selective pulses, in a 3-qubit register. 
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FIG. 2: Circuits for controlled gates U among two nuclear 
spins (A) and a nuclear and the electronic spin (B). The gates 
A, B, C are defined such that U = e ia AZBZC and ABC = 1, 
where Z is a 7v rotation around z 26]. $ a is a phase gate: 
|0> <0| + |l)(l|e iQ/2 and the gate a indicates e ia/2 t. 



fine couplings and Tidip the nuclear dipolar interaction, 
whose strenght can be enhanced by the hyperfine inter- 
action [ijj]. When the electronic spin is in the m s =l 
state, nearby nuclei are static (since distinct hypcrfinc 
couplings make nuclear flip-flops energetically unfavor- 
able in the so-called frozen core [3]) and give rise to 
a quasi-static field acting on the electronic spin. The 
other bath nuclei cause decoherence via spectral diffu- 
but their couplings, which determine the 



17, \U 



noise strength and correlation time, are orders of magni- 
tude lower than the interactions used to control the sys- 
tem. While in the m s =0 manifold all the nuclear spins 
precess at the same frequency, the effective resonance 
frequencies in the m s = l manifold, lo{, are given by the 
hyperfine interaction and the enhanced g-tensor 13, ll9( , 
yielding a wide range of values. Some of the nuclear spins 
in the frozen core can thus be used as qubits. When fix- 
ing the boundary between system and environment we 
have to consider not only frequency addressability but 
also achievable gate times, which need to be shorter than 
the decoherence time. 

Control is obtained with microwave (/iw) and radio- 
frequency (rf) fields. The most intuitive scheme, per- 
forming single-qubit gates with these fields and two-qubit 
gates by direct spin-spin couplings, is relatively slow, 
since rf transitions are weak. Another strategy, requir- 
ing only control of the electronic spin, has been proposed 
W ' 



Zeeman eigenstates induces nuclear spins rotations about 
two non parallel axes that generate any single-qubit gate. 
This strategy is not the most appropriate here, since ro- 
tations in the m s =0 manifold are slow f25|. We thus 
propose another scheme to achieve universal control on 
the register, using only two types of gates: 1) One-qubit 
gates on the electronic spin and 2) Controlled gates on 
each of the nuclear spins. The first gate can be simply ob- 
tained by a strong /iw pulse. The controlled gates are im- 
plemented with rf-pulses on resonance with the effective 
frequency of individual nuclear spins in the m s =l man- 
ifold, which are resolved due to the hyperfine coupling 
and distinct from the bath frequency [23[. Achievable 



rf power provides sufficiently fast rotations since the hy- 
perfine interaction enhances the nuclear Rabi frequency 
when m s =l 24] . Any other gate needed for universal 
control can be obtained combining these two gates. For 
example, it is possible to implement a single nuclear qubit 
rotation by repeating the controlled gate after applying a 
7r-pulse to the electronic spin. Controlled gates between 
two nuclei can be implemented by taking advantage of 
the stronger coupling to the electron as shown in Fig. 
[2jA). As long as the hyperfine coupling is several times 
larger than the nuclear coupling, the scheme avoiding any 
direct nuclear interaction is faster. Although selectively 
addressing ESR transitions is a direct way to perform a 
controlled rotation with the electronic spin as a target, 
this is inefficient as the number of nuclear spin increases. 
The circuit in Fig. [21(B) performs the desired operation 
with only the two proposed gates on a faster time scale. 

When working in the m s =l manifold, each nuclear 
spin qubit is quantized along a different direction and 
we cannot define a common rotating frame. The evolu- 
tion must be described in the laboratory frame while the 
control Hamiltonian is fixed in a given direction for all 
the nuclei (e.g. along the x-axis) . This yields a reduced rf 
Rabi frequency Q = f2 r f y/ cos (pi 2 cos6\ 2 + sini^i 2 (where 
ipx} define local quantization axes in the m s =l man- 
ifold and f2 r f is the hyperfme-enhanced Rabi frequency). 
The propagator for a pulse time t p and phase ip is 



21] : switching the electronic spin between its two Ui,(Cl T { t tp) 



_ -i[ut p -(\—tl>)]az/2 -intptTi/2 -t(A-V>)ff S /2 



where {crj, o-y , erj} are the Pauli matrices in the local 
frame and A is defined by tan (A) = tan ip±/ cos 9± . An 
arbitrary gate U = R^(j)Rx((3)Rs(a) can be obtained 
by combining Ul with an echo scheme (Fig. [3j) , which 
not only refocuses the extra free evolution due to the lab 
frame transformation, but also sets the gate duration to 
any desired clock time common to all registers. Fixing a 
clock time is advantageous to synchronize the operation 
of many registers. This yields a minimum clock time 
T>47rxMax n3 ^{n- 1 + (^")- 1 }. 

In order to refocus the fast electronic-spin dephasing 
given by the frozen core nuclear spins, we need to embed 
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FIG. 3: Rf pulse scheme for 1 nuclear spin gate in the m s =l 
manifold. With t p =/3/Cl and fixing ipi=a-\ and if)2=n-\, the 
time delays are £i = -|- — t p — — ^p- and t2=-j — t-x + ^p-- 



the control strategy described above in a dynamical de- 
coupling scheme |27j without loosing universal control, as 
explained in the following. The electron-bath Hamilto- 
nian is given by Eq. ([1]) , where the index j now runs over 
the spins in the bath. Neglecting for the moment the cou- 
plings among nuclei, we can solve for the evolution of the 
electronic spin under an echo sequence. By defining Uo 
and Ui the propagators in the and 1 manifold respec- 
tively and assuming that the nuclear spins are initially 
in the identity state (high temperature regime), we cal- 
culate the dynamics of the electronic spin, p e (t) = [1 + 

(\0)(1\ f ee (t)+h.c.)} Pe {0), where /„(*) - Tr ' 



J] [1 - 2 sin 2 (6{) sin 2 {u{t/2) sin 2 {to L t/2) 



is the func- 



tion describing the echo envelope experiments nam. 
Since in the m s —0 subspace all the spins have the same 
frequency, j ' ee (2mr / 'u> l) = and the electron comes back 
to the initial state. Nuclear spin-spin couplings lead to 
an imperfect echo revival and ultimately to decoherence 
via spectral diffusion 17, 0|. The energy-conserving 
flip-flops of remote nuclear spins modulate the hyperfme 
couplings, causing the effective field seen by the electron 
to vary in time. The field oscillations can be modeled 
by a classical stochastic process. The overall evolution 
of the electronic spin is therefore due to two processes 
that can be treated separately as they commute: the 
echo envelope calculated above and the decay due to a 
stochastic field, approximated by a cumulant expansion 
[29l | . For a Lorentzian noise with autocorrelation function 

G(t) — Goe /Tc wc obtain a spin-echo decay oc e 3t = 
for t <^ t c (or the motional narrowing regime for r c Ct). 

By using dynamical decoupling techniques [3(| and se- 
lecting a cycle time multiple of the bare larmor preces- 
sion period it is possible to extend the life-time of the 
electronic coherence. Figure [4] shows how to combine the 
electron modulation with the sequence implementing spin 
gates. The effectiveness of these techniques relies on the 
ability to modulate the evolution on a time scale shorter 
than the noise correlation time. The noise of the electron- 
nuclear spin system is particularly adapt to these decou- 
pling schemes. Consider for example the NV center. The 
measured electron dephasing T 2 * time is about 1/xs in nat- 
ural diamonds [l9| , as expected from MHz-strong hyper- 
fine couplings. The intrinsic decoherence time T2 can be 
orders of magnitude longer (T2 ^ 600yus). This reflects 
the existence of a frozen region, where the spin flip-flops 
are quenched. The radius of this frozen core is about 2.2 
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FIG. 4: rf and fiw pulse sequence to implement a 1 nuclear 
spin gate while reducing the effects of a slowly varying elec- 
tronic dephasing noise. The black bars indicate 7r-pulses, 
while the first rectangle indicate a general pulse around the 
x-axis. r = ? — % and t w = ^r^. 



nm [16[ and only spins with hyperfme coupling <2.5 kHz 
contribute to spectral diffusion. Both the inverse cor- 
relation time (given by the dipolar coupling among car- 
bon spins) and the noise rms (given by the coupling to 
the electron) are of order of few kHz. Dynamical decou- 
pling schemes must thus act on time scales of hundreds 
of fis. This in turns sets achievable constraints on the 
gate speed. 

The time of a conditional single nuclear spin rota- 
tion must be set so that the Rabi frequency f2 r f is much 
less than the frequency splitting between two neighbor- 
ing spins (in terms of frequencies). Suppose we want 
to control an n-spin register without exciting bath spins 
at the Larmor frequency. The minimum frequency split- 
ting between two nuclei in the m s =l manifold will be 
at best 5uj = for nuclear frequency equally spaced 
and u>m the maximum hyperfme- induced effective nuclear 
frequency. The nuclear frequency spread due to the hy- 
perfme interaction is then AE^ = ^^-ojm- We want 
the Rabi frequencies to satisfy the constraints: AE e 
n e > ^ujm > ^ > rf , where AE e = 2gp B B z is 
the gap to other Zeeman levels {m s =— 1 for the NVc) 
and fi e the fsw power. For a typical choice of 700G mag- 
netic field along the NV axis, we have AE e = 2GHz and 
ujl = 0.8MHz. Assuming u>m ~ 20MHz and n — 4 spins, 
we obtain the following parameter window (in units of 
MHz) 2000 > Q e > 25, fi rf <C 5. The gate clock time 
can be as short as a few [is, allowing hundreds of gates 
in the coherence time. 

Since the scheme presented is based on selective pulses, 
the most important (coherent) errors will be due to off- 
resonance effects. If the Rabi frequency is much smaller 
than the off-set from the transmitter frequency 5u>, the 
off-resonance spin will just experience a shift (Bloch- 
Siegert shift) of its resonance frequency, Ao; BS = Su> — 

>. This results in 



5uJ 2 dt >^_^ = _ 



an additional phase acquired during the gate time that 
needs to be refocused. Note that this error grows with 
the register size and constrains f2 r f. When reducing the 
Rabi frequency to achieve frequency selectivity, we have 
to pay closer attention to the rotating-wave approxima- 
tion and consider its first order correction, which pro- 
duces a shift of the on-resonance spin Aw RWA = f2 2 f /4wM- 
Other sources of errors are evolution of bystander nuclear 
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spins and couplings among spins. More complex active 



decoupling schemes [3JJ, |32|, |33[ can correct for these er- 
rors, allowing to use the m s ~l manifold as a memory. 

Advanced techniques like shaped pulses, with ampli- 
tude and phase ramping, composite pulses [34(, pulses 



optimized via optimal control theory or with numerical 
techniques [35|, 3f| 37 1 can provide better fidelity. The 
analytical model of control serves then as an initial guess 
for the numerical searches. Pulses found in this way cor- 
rect for the couplings among nuclei and are robust over a 
wide range of parameters (such as experimental errors or 
the noise associated with static fields). Table |T] shows the 
results of simulations in a fictitious NV system with 1-4 
nuclear spins and effective frequencies in the m s =l man- 
ifold ranging from 15 to 2MHz. We searched numerically 
via a conjugate gradient algorithm for a control sequence 
performing a desired unitary evolution, by varying the 
amplitude and phase of the fj.w and rf fields. We then 
simulated the control sequence in the presence of noise, 
with contributions from a large, static field and a smaller 
fluctuating one. The projected fidelities in the absence 
of experimental errors are very high, a sign that the fast 
modulation of the electron evolution effectively decou- 
ples it from the bath. The pulse robustness with respect 
to the noise is slightly degraded as the number of spins 
increases: the noise induces a spread of the electron res- 
onance frequency, and it becomes more difficult to find 
a sequence of control parameters that drives the desired 
evolution in this larger Hilbert space for any of the possi- 
ble electronic frequency. The fidelity degradation is how- 
ever modest, and can be partially corrected by increasing 
the control field intensity. Furthermore, combining these 
pulses in a dynamical decoupling scheme would provide 
an efficient way to coherently control the registers. 





1 spin 


2 spins 


3 spins 


4 spins 


F (ideal) 


0.9999 


0.9999 


0.9992 


0.9995 


F (noise) 


0.9994 


0.9995 


0.9975 


0.9925 


time 


5.0 /xs 


5.5 /xs 


6.0 /xs 


8.5 fis 



TABLE I: Simulated fidelities (F = |Tr [UtU] /2 N \ 2 ) at the 
optimal gate time in the presence of noise. The gate is a n/2 
rotation about the local x-axis for qubit 1 in a register of 1- 
4 nuclear qubits. The noise parameters are T 2 * = 1.5/xs and 
T2 — 250/is; the maximum Rabi frequencies are 2ir x 10MHz 
and 2ir x 20kHz for the electronic and nuclear spins respec- 
tively. As expected, the optimal gate time increases with 
the register size, reflecting both the more complex control re- 
quired in a larger Hilbert space and the weaker hyperfine cou- 
plings to more distant spins. Similar fidelities were obtained 
for a CNOT gate between spin 1 and 2. 

The size of the register is eventually limited by the 
number of available nuclear spins with a hyperfine cou- 
pling strong enough to be separated from the bath. From 
experiments and ab-initio calculation [38| we expect hy- 
perfine couplings of ~ 130MHz in the first shell, and then 



a number of about 50 possible nuclear sites with hyperfine 
values from 15MHz to 1MHz. Even if the concentration 
of C-13 is increased (and the Nitrogen nuclear spin is 
used) the size of the register will be limited to about 10 
spins. Nevertheless such registers would be very useful 
for memory storage and error correction purposes. 

In conclusion, we have presented a general approach 
to the control of a small quantum system comprising an 
electronic spin and few nuclear spins in the surround- 
ing spin bath. We have shown that several of the bath 
spins can be isolated and effectively controlled, yielding 
a few-qubit register. These registers can be employed 
in many proposals for distributed quantum computa- 
tion and communication, where coupling among registers 
could be provided either via photon entanglement [HI 
or by a movable reading tip [391 ]. Our control methods 
enable algorithms needed for error correction and entan- 
glement purification, while the nuclear spins provide a 
long-time memory in the m s — 1 manifold, via active 
refocusing, and the electron dephasing is kept under con- 
trol by dynamical decoupling methods. We thus develop 
a modular control scheme, scalable to many registers and 
applicable to many physical implementations. 
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